Novel Mechanisms for repulsive Casimir forces 
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We present two novel models for repulsive Casimir interaction between positive perturbations. 
One example relies on non locality of the dielectric response and one relies on interference between 
(attractive) modes. Such examples are impossible to achieve in Id massless theories, as they are 
prohibited by a generalization of the result of Lambrecht, Jaekel and Reynaid [lj to non local 
dielectrics and to negative perturbations. 
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I. INTRODUCTION 

One of the intriguing aspects of the Casimir effect is 
it's subtle dependence on geometry and type of bound- 
aries. A lot of effort and progress has been made recently 
in quantifying the force in various situations, mostly with 
experimental applications in mind utilizing scattering ap- 
proximations as well as numerical ones. From both prac- 
tical and theoretical point of view, reversing the sign of 
the force from attractive to repulsive is of particular inter- 
est. This requires considering non-symmetric situations 
to avoid the theorem 0, [j| ensuring attraction between 
opposing bodies having the same properties. The depen- 
dence of the sign on asymmetric boundary conditions was 
studied in many works, most recently, in the context of 
critical Casimir forces [6j where the direction of the force 
for various boundary conditions was explored for inter- 
acting theories under a renormalization group flow. 

The intuition regarding the sign of Casimir forces may 
be most easily understood from second order perturba- 
tion theory. This includes pair-wise summation for di- 
electric bodies, or the two scattering approximation of 
[Til ] . The result of this approach is attraction for any two 
dielectrics in vacuum. On the other hand, perfectly mag- 
netically polarizable and perfectly electrically polarizable 
atoms will repel each other [7j, with the consequence of 
repulsion between a perfectly conducting and perfectly 
permeable materials 

SUES- 

For many purposes it is convenient to "soften" the 
boundary conditions, by considering the interaction of 
the free field with the material, as a perturbation on the 
free field wave equation. Consider for example, two di- 
electrics A and B immersed in a medium M, of dielectric 
constant tu- The wave equation in this case is written 
as 



affect the mutual force) 



V x V x A + uj 2 e(x)A = 



(1) 



where uj 2 e(x) = io 2 (-M + Va + Vb is a sum of the constant 
background term w 2 £m, and two perturbation terms 
Va(x) = u 2 e A {x) - LL> 2 e M ,V B (x) = uj 2 e B (x) - uj 2 e M rep- 
resenting the response of the two bodies. 

Second order perturbation theory for the energy of the 
system gives (ignoring self energy terms which do not 
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Tr(V A G V B G ) 



(2) 



More explicitly this may be written as the integral (see 
,1313) as 



e.g. [12 



E& cx ( 
-f^f A dxf B dyV A (x, iu)V B {y, iu>) (G 0ij (x - y)) 

where G (x,x% = (x| 7x7x 1 +u2tM |x%- is the Euclid- 
ian (Wick rotated) propagator in the medium M. Since 
Go(x,x')y is real, Eq. © shows that the sign of the in- 
teraction energy depends only on the signs of the pertur- 
bations Va,V b as sgn[E^] = —sgn[VA]sgn[V B ]. Since 
E^ vanish for infinitely separated bodies, an immediate 
conclusion is that at least at large distances, perturba- 
tions of similar signs attract while perturbations of op- 
posite sign repel. It seems natural to conjecture this to 
hold also at arbitrary distances. This however does not 
follow from the above argument, since the energy need 
not be a monotonic function of the distance. 

Eq. ((3]) is an easy way of understanding the well known 
situation where the Hamaker constant is negative, when 
considering materials A, B such that e A < £m < (b ■ 

Using the more abstract expression Eq. ^ and the 
fact that Go is a positive operator allows generalizing 
the above arguments to the case where the perturba- 
tions V Al V B are given by arbitrary operators of defi- 
nite sign. (An hermitian operator V is said to be pos- 
itive/negative iff all its nonzero eigenvalues are posi- 
tive/negative. This is equivalent to requiring (?/>|V r |'(/>) > 
for all states 

One example of such generalized perturbation is to as- 
sume nonlocal dielectric functions. As another example 
imagine the bodies A, B to have nontrivial magnetic per- 
meability. In this case Eq. (JTJ) would have to include the 
extra term V x (^jy — 1)V X A which may be interpreted 
as a perturbation VA with V a differential operator. For 
fi > 1 one may show that this differential operator is 
negative. Thus the result @, H, @, [HJ. that a dielec- 
tric wall repel from a permeable one is consistent with 
the mnemonic that "positive" perturbations attract one 
another while positive perturbation will repel a negative 
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All known examples of repulsive Casimir forces can be 
traced to such a mechanism. One may thus expect that 
the sign of the Casimir force is the product of the signs 
of the perturbations, whenever these perturbations have 
a definite sign. 

There are two general statements regarding the sign 
of the force at all distances: In Id, Lambrecht, Jaekel 
and Reynaud [l[ have shown that dielectrics attract un- 
der very general conditions. We shall show below that 
this statement also hold for general positive perturba- 
tions in Id. Unfortunately, this argument does not ex- 
tend to higher dimensions. The only available result to 
date is that under conditions of mirror image symmetry, 
the force is always attractive 0, [|| . 

A natural question that rises is: is attraction always 
guaranteed between a pair of positive perturbations or a 
pair of negative perturbations? If not, what other mech- 
anism can be responsible for repulsion ? 

In this paper we give a partial answer to the question 
raised above: We show that this intuition is, in fact, cor- 
rect in ID, i.e. any two perturbations (in the general 
sense) having the same sign attract, by which we gener- 
alize the result P, . On the other hand, the situation 
in higher dimensions is not as simple: In fact, we find 
examples of positive perturbations which yield repulsion. 
These examples, do not rely on the relative sign of the 
perturbation as previous examples do, and thus provide 
new mechanisms for repulsion. While the examples we 
present are "toy" models, it is important to stress that 
both examples may be thought of as extreme limits of 
physical perturbations. 

A convenient starting point for our analysis is the 
TGTG representation for the energy Q: We consider 
the interaction between two bodies A and B which are 
defined through their dielectric susceptibilities xa = 
tA — 1, xb — £b — 1 respectively or any generalized per- 
turbation (i.e. xa can be an differential or integral oper- 
ator). The regularized Casimir energy is given by: 



E c {a) 



f °°^logdet(l 



TaGqTbGo) 



(4) 



Here Go is the free (Hclmholtz) Green's function, given 
by G (x, x') = (x| _y} +u -2 |x'), for a scalar field (or 
G (x,x% = (x| VxV 1 x+ ^ |x%- when considering the 
electromagnetic field). The operators T^,Tb are a Wick 
rotated version of the the transition operators appearing 
in the Lippmann-Schwinger equation. They are given by 



A.B. 



In general, the determinant formula for the Casimir en- 
ergy is hard to analyze because it involves the determi- 
nant of an infinite dimensional operator. However, the 
fact that the TGTG determinant in (jlj is well defined , 
means that in practice only a finite dimensional subspace 
gives significant contribution to the force. Such a finite 



dimensional reduction allows arbitrary good numerical 
evaluation of the force. This suggests that the question 
of the possible sign of the force may also be addressed 
within the context of finite dimensional operators. In 
any case, to get the correct sign of the force, one should 
be careful when using infinite dimensional tricks such as 
the "identity" = -i < 0. 

In this paper we concentrate on two situations where 
the determinant in (TJJ reduces to a strictly finite dimen- 
sional one. The first is the case of the so called separa- 
ble potentials. These are (typically non-local) potentials 
which by construction are given by finite rank operators 
and therefore interact only with a finite dimensional sub- 
space of states. The second is the case of one dimensional 
systems. At each given energy to a one dimensional field 
(j> has only two modes: the left and the right mover state 
|Z),|i2). (A multicomponent <fi = (0i, fa, ■■■fa) would 
have 2n modes). The determinant (|4|) then becomes an 
ordinary finite dimensional one. 

Below, we show that both of these models allow simple 
examples where the interaction energy between a pair of 
local positive perturbations may be computed exactly. 
We then show how the perturbations can be chosen in 
a way to yield a non-monotonic energy dependence on 
distance, implying that the force can reverse its direction, 
even when both constraints are positive. 



II. ATTRACTION IN ID FOR GENERAL 
POSITIVE PERTURBATIONS 



The purpose of this section is to establish that in ID, 
any two positive perturbations attract each other. In one 
dimensional systems the TGTG formula for the Casimir 
interaction between objects A, B situated a distance a 
apart reduces to well known relation (l6| : 



E 



2^ 



log(l-e- 2a ^ A Mr B M), 



(5) 



where rj^iui), rs(iij) are reflection coefficients for plane 
waves scattering on the bodies A, B, evaluated on the 
imaginary frequency axis. It was shown in [l| that for 
a material whose dielectric function is local, r(iu>) < 
0, Voj € R, which has the immediate consequence that in 
one dimension any two local dielectrics attract. Here we 
show that this result in fact applies to any kind of positive 
perturbation, including the important case of non-local 
dielectric functions. 

To do so, we first write the reflection coefficient from 
a Id potential as the matrix element 

r(u) = -^(R\T\L) = / dxdx'{e' Mx yT(x,x' : uj)e luJX \6) 

Thus after Wick rotation we have r(iui) — 
J dxdx'e~ u ' x T(x,x' ,iui)e~ UJX . (The integral converge 
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since we consider only potentials of compact support 
and hence also T is of compact support.). It has been 
shown [l3l |, that a positive potential V > 0, implies 
that at imaginary frequencies the Lippmann-Schwinger 
operator is positive T(ilu) > [21] . Note that r(iuj) is 
of the general form — (ip\T(iuj)\ip} so that the positivity 
of T(iuj) implies that r(iuj) < 0. 

Similar arguments may be applied to a Id multi- 
component field <fi — ^2, ■■■4>n), assuming all compo- 
nents are massless (or have exactly the same mass) . The 
reflection coefficients ta,b in this case turn into n x n 
matrices which by the argument given above are strictly 
negative matrices. The energy may be expressed as 

E = J ^ log det (1 - e-^f A (ito)r B (tw)) , (7) 

where the determinant is an ordinary n x n one. Since 
spec(r A r B ) = spec(v '-r A {-r B )^ '-r A ) = {Ai,A 2 ,..A n } 
is positive, the energy can be written as 



E = / ^£log(l-e- 2 -A fe (u,)) 

k=0 



(8) 



where < Ai, A2, ..An < 1 are a-independent. (The in- 
equality < A < 1 has been shown in [l3j]) We immedi- 
ately see that: 



F, 



c 



d a E c 



/duj s—^ 2ue 



-2au 



Afc(w) 



- 2a "A fc M) 



<0 (9) 



Alternatively, one could arrive at the same conclu- 
sion directly from ([4]) by noting that in the Id scatter- 
ing basis Go appearing in the expression TGTG is es- 
sentially the c-number e~ auJ . Since positivity of Va,V b 
guarantee that Ta,T b are also positive operators, as 
shown in Q it is enough to have required property 



of TGTG 



-2ka 



T A T B . However the spec(T A T B ) 



spec(v7A^B \Ta) is positive, therefore the energy can 
be written as 



J l n 



A fc H) 



(10) 



where < Ai, A 2 



fe=0 

.A„ < 1 are a-independent. 



III. INTERFERENCE INDUCED CASIMIR 
REPULSION FOR FIELDS WITH DIFFERENT 
DISPERSION RELATIONS 

Formally, a field in a space of arbitrary dimension may 
be thought of as a field in one dimension by consider- 
ing the transverse momentum k± as if it was an internal 
continuous index. Thus one may ask why doesn't the Id 
result ((9]) extend to higher dimensions. 

The answer has nothing to do with the fact that there 
are infinitely many k± modes. Rather, it is related to 



the fact that in strict Id we have (implicitly) assumed 
that in free space (i.e. outside the scatterers) all the 
component of <f> — (cf>x, <p2, •■■4> n ) share the same (zero) 
potential. This is in contrast to the higher dimensional 
problem where the term k\ in the lagrangian serves to 
distinguish between the different scattering channels even 
in empty space and makes the fields in the Id picture have 
different masses. As a result, the "Id momentum" fen is 
no longer conserved, which from the Id point of view is 
a source of great extra complications. 




FIG. 1: Bodies A and B in a confined space, transverse modes 
are quantized. In this situation, only low transverse modes 
contribute. The problem becomes equivalent to a finite num- 
ber of Id fields with different masses. 

A simple toy model corresponding to the above pic- 
ture may be constructed by taking <f) = (4>ii4>2, ■■■<Pn) 
to be Id fields of masses m 1 ,m 2 ,...m„. Let us denote 
fcj = \J uj 1 + mf and let K be the diagonal matrix hav- 
ing ki as its eigenvalues. The (imaginary frequency) free 
propagator may then be expressed as Go = 2lf e • 

Consider enforcing on cj) the Dirichlet type boundary 
conditions: J2 a i4>i\x=a = 0, Y^,Pi<l>i\x=b = 0. Here a = 
(ai, 0L2, ■■■■ot n ) and (3 = (j3\, fa, ■■■■f3 n ) are some constant 
vectors. This may be obtained by assuming the (positive) 
potential V = Va + V B 

Va = Act <S> a^8(x — x a ), V B — A/3 ® ^5(x — Xb) 

with A — > 00. One then finds the reflection coefficients 
from the two potentials: 



r B (iuj) = r B (iuj) 
where we used the notation 



|d| 2 +2w/A' 
|/3| 2 +2w/A 



(11) 



o 



'0. 



and similarly [3 = 
We note, parenthetically, that a more general dispersion 
require — > i.e. the group velocity is the extra 

factor. Calculating the determinant one obtains 



det(l 



rAe -K\x a -x b \f Be -K\x a 



-x b \\ 



) = 



(12) 



1 - 



1/5) 



|fi|.|/3| 



It is quite easy to see that this expression need not be 
monotonic as a function of the distance x — \x a — Xb\- 
To check the behavior of the total (w-integrated) energy 
is more complicated. However by choosing special values 
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for the parameters one may check numerically that the 
resultant energy need not be monotonic. The following 
graph, Fig. shows the interaction energy for the case 
ri = 2 with 7712/mi = 02/01 = —(P2/P1) = 5 • 
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FIG. 2: Casimir energy as a function of distance (measured 
in units where mi = 1), obtained from integration of Eq. (|12p 
over frequencies. Grey area marks the local energy minimum. 

We stress that each of the field modes <f>\ and <p2 on 
each own would produce an attractive force if the other 
was somehow turned off (e.g. by giving it very large 
mass). The repulsion seen in Fig. [5] may therefore be 
understood as an interference between the modes (or, 
more precisely, from the dependence of the interference 
term on the separation). As a side remark, we mention 
that the mathematics involved here is somewhat reminis- 
cent to that of the Glashow- Iliopoulous-Maiani (GIM) 
mechanism 17|, although actually applying it to quark 
mediated weak interactions would be, of course, com- 
pletely unrealistic. 

For a realistic way to obtain such a model, consider e.g. 
objects situated inside a circular cylinder of radius R as 
in Fig. [T] Standard mode expansion of a field leaving 
in this cylinder would reduce it to a series of Id fields of 
masses m n = Cn/R with Cn,m zeroes of Bessel functions 
J m and of each derivative for TM and TE modes re- 
spectively. Standard boundary conditions on the bodies 
would demand specific (shape dependent) combinations 
of the Id fields to vanish. 

The fact that for very large and for very small dis- 
tances the force attracts is a general feature not special 
to the specific model. At very large distances x ^> 1 
only the channels of lowest mass contribute to the force. 
Having only this single mass, the above no-go theorem 
applies and we obtain attraction. At very small distances 
x <C 1 the integrand decays very slowly with uj so that 
most of the contribution to the force comes from very 
large frequencies w ~ A. At such frequencies all the fcj's 

*> u) Vi and 



become practically equal: fcj = yj m| + a, 
our theorem again imply attraction. 

(E.g. in the above model E ~ -^Polylog[2, 
At intermediate distances the force may however change 
sign as the above example demonstrates. 

It should be remarked however that this reasoning re- 
lies on the assumption of having only finitely many chan- 
nels. In an infinite channel system it is conceivable that 
as x becomes smaller more and more new channels (of 
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FIG. 3: An effective interaction with a separable potential can 
happen when a free field is perturbed by tunnelling amplitude 
into traps, A and B in which the field may have bound states. 



higher m) become relevant, making the above argument 
false. 

The short distance attraction argument may also be- 
come false in case the reflection coefficients decay very 
rapidly as w — > 00. Also we remark that using fine tuning 
one might be able to construct a model where the leading 
short distance attractive term vanishes, in which case the 
resultant force can remain repulsive down to x = 0. For 
example: Suppose one constructs a model where the o's 
and /3's depend on to in such a way that a, (3 are constants 
and take (n = 2) 02/01 = —P2/P1 = = 0,TO2 7^ 

then the attractive term which usually dominates at 
x — > vanishes and the next to leading term result in 
a force which is attractive down to x = 0. However any 
slight change of the parameters will make the usual lead- 
ing term re-appear and thus turn the force into an at- 
tractive one (at x = 0). 



IV. NON-LOCALITY: REPULSION BETWEEN 
SEPARABLE POTENTIALS 

Our next repulsive example shows that for d > 1 non- 
locality of dielectric functions, may also, in principle, re- 
sult in repulsive behavior. The effect of spatial dispersion 
on Casimir forces between metals has been of interest for 
some time 16|, . It is particularly important in metals 
where long range density correlations in the medium can 
be present and substantially change the scaling of the 
force [H. 

Our example here is very different from the situa- 
tions considered in the references above. Our toy model 
is based on so called "separable potentials". Such po- 
tentials arise in variety of situations in physics and in 
mathematics, and were first introduced for Casimir type 
problems by Jaffe and Williamson JlS]. They corre- 
spond to "rank 1" perturbations, and can be written as 
V = I/) (/| for some function /, or in x space notation 
as V(x,x') = f(x)f*(x'). Such separable potentials de- 
scribe the interactions involved in Feschbach resonances, 
and play a major role in the properties of fermi conden- 
sates. Note, also, that such V is a positive operator by 
construction since (ip\V\ip} — \(^p\f}\ 2 > for any if). 

The situation we consider is of a scalar field, 
which can tunnel into traps where it has a bound 
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state Fig. [3l Thus, in our case the perturba- 
tions are traps. Integrating out degrees of freedom in 
these traps will generate a perturbation of the form 
J dxdx' 4>(x)t(x)Gf > ound(x , x' )t* (x')4>(x r ) for the field out- 
side the trap, where Gbound is the green's function of the 
trap, and t(x) is the tunnelling clement. 

Close to resonance, Ground may be approximated as: 

, fb(x)fb(x') 



G 



E-E h 



where ft is the wave function of a 



bound state, and when E is close to E^. Thus, the effec- 
tive interaction with the trap is through the separable po- 
tential V = f(x)f*(x'), where f(x) = J \ E l Eb \ t{x)f b {x). 

For such potentials, one may readily calculate the T 
operators, which turn out to be also of rank one, and 
so the interaction energy is easily calculated (One may 
also look at sums of separable potentials, corresponding 
to finite rank perturbations, which can be computed in 
the same way). 

In case the free field interacts with the perturbation 

era tut 

oc !/)(/!• We have: 



T 



l + VGo 



V 



l + </|C?ol/> 



i+WGoW 

l/X/l 



(13) 



V. SUMMARY 

To summarize, we gave a general rule for the sign of the 
Casimir force at large distances and considered mecha- 
nisms allowing the force to reverse its direction at shorter 
distance. We have shown that such behavior is possible 
when an interaction is introduced between fields with dif- 
ferent dispersion relations. Such a situation may mimic 
the interaction between bodies confined in an infinite 
cylinder which serves to quantize the transverse modes, 
rendering a quasi Id situation. Secondly, we have shown 
that non-locality of a dielectric function may, in an ex- 
treme limit, result in repulsive behavior. This type of 
non-locality, typical of separable potentials, may arise in 
situations where the bodies represent "traps" into which 
the field can tunnel, and interact with a bound state. 
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And so, 



log(l 



logdet(l - T A G T B G ) = (14) 
1 :\(fB\G \f A )\ 2 ) 



1+(/a|G |/a) 1+(/b|<?o|/b 



thus, since the (positive) terms 1+{f J GolfA) 1+{fB l G „\f B ) 
do not depend on the distance, to find the direction of 
the force it is enough to consider |(/b|Go|/a)|- 

As an example, we consider the following functions 

fA = (1 + ^+i)[ai5{x + a) + a 2 S{x + a + 0.1)] (15) 
Sb = (1 + J+i)[Pi8{x ~a) + P 2 S(x - a - 1)]. 



FIG. 4: Casimir energy as a function of distance a between 
separable potentials Va = |/a)(/a| and Va = |/s)(/s| for 
/a, /b denned in (Tl5)) . 
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